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Abstract. We compare population synthesis results for inner and outer magnetosphere emission 
i— i models with the various characteristics measured in the first LAT pulsar catalogue for both the 

radio-loud and radio- weak or radio-quiet y-ray pulsars. We show that all models fail to reproduce 
the observations: for each model there is a lack of luminous and energetic objects that suggest a 
non dipolar magnetic field structure or spin-down evolution. The large dispersion that we find in 
the simulated gamma-ray luminosity versus spin-down power relation does not allow to use the 
present trend seen in the Fermi data to distinguish among models. For each model and each Fermi 
detected pulsar, we have generated light curves as a function of obliquity and inclination angles. The 
theoretical curves were fitted to the observed one, using a maximum-likelihood approach, to derive 
the best-fit orientations and to compare how well each model can reproduce the data. Including the 
radio light-curve gives an additional key constraint to restrict the orientation space. 
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RADIO PULSAR POPULATION SYNTHESIS AND EMISSION 
°i MODELS 
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We synthesised a population of 10 1 neutron stars (NSs). The stars were born on the 

galactic plane with a surface density distribution obtained from the location of the HII 

clouds in the Milky Way and spread in the Galaxy by assuming a maxwellian supernova 

kick velocity distribution with mean of 400 kms -1 . Each star has been spun down as 

a rotating dipolar magnetic field and by assuming a constant exponential rate magnetic 

field decay on a timescale of 2.8 Myr. For the details and references on the NS sample 

evolution see the Gonthier et al. proceeding of this same conference. 
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Gamma-ray and radio emission models 

For each NS of the simulated sample we have evaluated y-ray luminosity and emission 
pattern according to 4 different emission models: Polar Cap, PC (low magnetosphere, 
two pole caustic) [1], Slot Gap, SG (intermediate-high magnetosphere, two pole caustic) 



[2], Outer Gap (OG) [3] and One Pole Caustic (OPC) [4] (both high magnetosphere, one 
pole caustic). The radio emission from each NS of the sample has been simulated by 
assuming an emission cone composed by a core and cone components [5]. 

Gamma & radio visibility 

For each assumed emission model, the visible component of the whole pulsar popu- 
lation has been obtained by applying a one year LAT pulsar visibility map. Each y-ray 
visible subsample has then been scaled to the number of radio pulsars in 10 radio sur- 
veys of the ATNF catalogue. This gives the number of simulated visible pulsars that 
would have been observed by the LAT telescope on the FERMI satellite after one year 
of observation, by assuming each emission model. 

FITTING THE PULSAR EMISSION WITH THE IMPLEMENTED 

MODELS 

By using the emission pattern of the implemented emission model, we performed an 
estimate of the magnetic obliquity a (angle between magnetic and spin axis) and line 
of sight £ (angle between spin axis and observer direction) for a sample of 22 radio 
loud (RL) LAT pulsars. The estimate has been implemented by fitting jointly the y- 
ray and radio light curves. A second (a; Q estimate has been performed by fitting just 
the y-ray LAT pulsar light curves with the high-energy model emission patterns. The 
comparison between the (a; Q solution space of the joint radio+gamma and gamma fit 
will be discussed in the results section. In both the fit methods and for each LAT pulsar, 
we have found a fit solution for each one of the assumed emission models. For each 
fitted pulsar, the model fit solutions (PC, SG, OG or OPC) characterised by the highest 
significance (highest likelihood value), is called best fit solution. 

RESULTS AND DISCUSSION 

Comparison LAT observations model predictions. In figure 1 are shown the com- 
parisons between the LAT observations and simulations, for luminosity vs. spin-down 
luminosity (left) and spin-down luminosity distribution (right). The PC, SG, and OPC 
model predictions well describe the expected trend L 7 oc E 05 while the OG model shows 
a too high dispersion. All the models show a lack of luminous objects and so fail in de- 
scribing the observed population. The gamma luminosity comparison cannot be used 
to discriminate which model best describe the observations. Figures 1 right shows that 
there is a lack of energetic objects in the model predictions. No one of the proposed 
models is able to explain the high number of high E pulsars observed by the LAT. The 
lack is obviously non model dependent and is not due to the choice of intrinsic param- 
eters of the stars, like radius, mass or moment of inertia. Several configuration of these 
intrinsic parameters were tried and all of them led to the same lack of energetic objects. 
The current mass, radius, and moment of inertia configuration has been chosen to boost 
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FIGURE 1. left: for each model is shown the comparison between L r = f(E) simulated (colored) and 
observation (black), right: E distribution comparisons between data (black) and simulations In both the 
plots, pink and green respectively refer to RL & RQ simulated pulsars. 



the power and minimise the observed lack. We identify two most probable explanations 
to the lack of energetic objects. A first one is based on a different magnetic field struc- 
ture (non dipolar) and so spin-down luminosity evolution for very energetic objects. A 
second one hypothesises a magnetic obliquity a evolution during the first stages of a 
pulsar life. The different shape between simulated and observed E distributions (fig. 1 
right) could support the non dipolar spin-down luminosity evolution scenario for the 
most energetic pulsars. The fact that the best agreement between data and simulations is 
obtained for the OPC model and since it differs from the OG just for the dependency of 
the luminosity from a (OG: L 7 = f(E, a); OPC: L r = f(E)\ could be one element in 
favor of the the magnetic obliquity evolution scenario and so, a possible a evolution for 
energetic objects should be taken into account. 




FIGURE 2. (oc;Q solutions plane for the joint radio plus gamma fit results (left) and just gamma results 
(right). In the right plot, triangles and square respectively refer to radio quite (RQ) and radio loud (RL) 
pulsars. In both the plots, pink, red, blue and black refer to PC, SG, OG, and OPC model solutions, and 
empty and filled markers respectively refers to all the solutions and to the best fit ones. 



Fit of the LAT pulsar profiles. In figure 2 and 3 are shown two of the main and 
suggestive results of the joint fit estimate of the LAT pulsars orientation. Figure 2 shows 
the comparison between the (oc;Q solutions plane for the joint radio plus gamma and 



just gamma fit. The gamma fit solution plane (figure 2 right) contains RL solutions 
in a region of the plane where is not possible to have radio emission. As it is shown 
in figure 2 left panel, to be able to explain both the radio and gamma emission, the 
RL pulsar, can have (oc\Q solution just along the radio diagonal of the a — £ plane 
(a;£ = 0;0 -> a;£ = 90; 90). This is justified by the fact that to be able to observe 
simultaneously radio and gamma emission, the angular distance between the observer 
line of sight £ and the magnetic axis can not be bigger of the angular radius of the radio 
emission beam that, for young pulsars, assume the values of ~ 20 — 30 degrees. 
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FIGURE 3. High energy cutoff as a function of the accelerator gap width. Green and black lines 
respectively refer to the best fit power law for all the solutions and for the highest significance ones. In both 
the plot, empty and filled markers respectively refers to all the solutions and to the highest significance 
ones. 

In figure 3 is shown the relation between the high energy cutoff and the gap width we 
found, for each model, in the fit results and not in the population study. This relation is 
particularly important because connect an intrinsic pulsar parameter like the with of the 
accelerator gap and the observed high energy cutoff. From the SG models, we obtained 
that E a* oc Wg°p 5 . 

Conclusions. The comparison between luminosity evolution and E distribution of 
simulated and observed samples show a lack of energetic objects in the model pre- 
dictions. As a possible cause, different spin-down luminosity evolution/magnetic field 
structures and/or a evolution have to be more investigated. The comparison between 
the (a; Q plane solutions of the joint fit and gamma fit method suggest that a multi- 
wavelength pulsar light curve fit is the only way to obtain reliable (a; Q estimates. The 
energy cutoff-gap width relation found offers the possibility for a future improvement 
of the population synthesis. 
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